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ABSTRACT 

The  propagation  of  electromagnetic  radiation  through  the  atmosphere  is  a  crucial  aspect 
of  laser  target  acquisition  and  surveillance  systems  and  is  vital  to  the  effective  implementation 
of  some  Theater  Missile  Defense  systems.  Atmospheric  turbulence  degrades  the  image  or  laser 
beam  quality  along  an  optical  path.  During  the  past  decade,  the  U.S.  Air  Force's  Geophysics 
Directorate  of  Phillips  Laboratory  collected  high  speed  differential  temperature  measurements 
of  the  atmospheric  temperature  structure  parameter,  Cj.2,  and  the  related  index  of  refraction 
structure  parameter,  Cn2.  The  stratospheric  results  show  a  1-2  order  of  magnitude  increase  in 
day  turbulence  values  compared  to  night.  Resolving  whether  these  results  were  real  or  an 
artifact  of  solar  contamination  is  a  critical  Theater  Missile  Defense  issue. 

This  thesis  analyzed  the  thermosonde  data  from  an  experimental  program  conducted  by 
the  Geophysics  Directorate  in  December  1990  and  found  strong  evidence  of  solar  induced 
artifacts  in  the  daytime  thermal  probe  data.  In  addition,  this  thesis  performed  a  theoretical 
analysis  of  the  thermal  response  versus  altitude  of  fine  wire  probes  being  used  in  a  new 
thermosonde  system  under  development  at  the  Naval  Postgraduate  School.  Experimental  wind 
tunnel  measurements  were  conducted  to  validate  the  analytical  predictions. 
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I.    INTRODUCTION 

The  propagation  of  electromagnetic  radiation  through  the  atmosphere  is  a 
crucial  aspect  of  many  laser  target  acquisition  and  surveillance  systems  and  is 
vital  to  the  effective  implementation  of  the  Theater  Missile  Defense  System  (TMD). 
Turbulence  degrades  the  quality  of  atmospheric  propagation  along  the  optical 
path  of  the  radiation.  Random  fluctuations  in  the  atmospheric  index  of  refraction 
produce  optical  turbulence  and  distort  the  phase  front  of  the  beam  [Ref.  1  :p.  6-6]. 
The  index  of  refraction  structure  parameter,  Cn2,  which  is  a  function  of 
atmospheric  wind  shear,  pressure,  and  temperature,  quantifies  optical  turbulence. 
Nearly  three  decades  of  measurements  show  that  Cn2  varies  both  horizontally, 
vertically,  and  with  time.  What  is  not  well-known  is  the  fine-scale  spatial 
distribution  and  the  intensity  of  turbulence  in  the  upper  troposphere  and  the 
stratosphere.  Atmospheric  temperature  fluctuations  characterized  by  the 
temperature  structure  parameter,  CT2,  produce  index  of  refraction  variations  by 
altering  the  air  density.  During  the  past  decade,  the  U.S.  Air  Force  Geophysics 
Directorate  of  Phillips  Laboratory  (PL/GP)  has  collected  measurements  of  this 
quantity  using  balloon-borne,  high-speed,  differential  temperature  thermosondes 
[Ref.  2:p.30].  Recently,  the  Naval  Postgraduate  School  (NPS)  has  used  both  fine 
thermistors  as  well  as  fine  wire  differential  thermocouples  to  measure  CT2  in  both 


day  and  night  conditions  [Ref.  3]  and  [Ref.  4].  The  GP  experiments  have  shown 
a  1-2  order  of  magnitude  increase  in  day  turbulence  values  compared  to  night. 
Attempts  to  explain  this  result  have  followed  two  lines  of  reasoning:  1)  either 
there  is  a  real  increase  in  optical  atmospheric  turbulence  in  the  day,  or  2)  the 
data  (collected  by  two  3.8  //m  platinum-coated  tungsten  wire  probes  for  GP  and 
12.4  /jm  chromel-constantan  thermocouples  for  the  NPS  Group)  are  unreliable 
in  the  day  due  to  systematic  effects  such  as  the  solar  heating  of  the  probes 
coupled  with  package  rotation,  pendulum  motion,  and  probe  vibration.  [Ref.  2:p. 
55] 

This  thesis  investigated  two  systematic  error  effects  in  balloon  thermosonde 
systems  used  by  NPS  and  GP.  This  was  accomplished  in  two  phases.  The  first 
was  to  look  for  solar  induced  artifacts  in  the  daytime  thermal  probe  data.  The 
second  phase  was  to  predict  then  measure  the  thermal  probe  response  time  of 
the  NPS  probes  as  a  function  of  altitude.  The  initial  phase  of  this  investigation 
consisted  of  a  statistical  evaluation  of  a  December  1990  experiment  conducted 
by  GP  at  Holloman  Air  Force  Base,  New  Mexico,  that  focused  on  the  cross- 
correlation  of  measurements  taken  from  adjacent  pairs  of  probes  under  day  and 
night  conditions.  The  second  phase  was  an  analytical  development  and 
measurement  of  the  probe  response  time  versus  altitude,  focusing  on 
thermodynamic  and  fluid  dynamic  factors.  These  calculations  were  compared 
with  wind  tunnel  measurements  to  validate  and  improve  current  probe  response 
models.     The  underlying  motivation  of  this  investigation  was  to  investigate 


systematic  error  sources  in  critical  atmospheric  turbulence  measurements  needed 
for  effective  deployment  of  theater  missile  defense  systems. 


II.   BACKGROUND 


A.      THEORY 

Lasers  emit  radiation  in  a  highly  directional,  collimated  beam  with  a  very  low 
angle  of  divergence.  The  beam  divergence  angle  9  can  be  "diffraction  limited," 
i.e.,  limited  by  the  wave  characteristics  of  light.  Under  vacuum  propagation 
conditions  the  divergence  angle  6  is 

e~^,  (i) 

where  A  is  the  light  wavelength  and  D  is  the  diameter  of  the  aperture.  [Ref.  5:p. 
21]  Atmospheric  turbulence  degrades  the  spatial  coherence  of  a  beam  of 
electromagnetic  radiation.  The  Navier-Stokes  equations  govern  the  fine-scale 
random  fluctuations  of  turbulent  flow  [Ref.  6:p.  262].  At  high  Reynolds  numbers, 
Re~105,  the  Navier-Stokes  equations  do  not  have  unique  solutions. 
Consequently  we  must  use  statistical  models  based  on  semi-empirical, 
dimensional  analysis  and  physical  reasoning.  These  focus  on  the  mean  flow 
properties  of  the  fluid  and  on  second  order  statistics  of  the  fluctuations 
(correlations,  spectra,  and  structure  functions)  [Ref.  7:  p.  311].  A  collimated  laser 
beam  initially  propagates  with  a  uniform  wave  front  until  it  encounters  a  turbulent 


medium  (the  atmosphere).  Atmospheric  temperature  inhomogeneities  cause 
random  fluctuations  in  the  atmospheric  index  of  refraction,  n,  [Ref.  1  :p.  6-6]  which 
distort  the  phase  of  the  wave.  The  result  of  this  distortion  on  an  initially 
collimated  laser  beam  is  to  induce  centroid  wander,  broadening,  and  break-up 
(in  strong  turbulence)  which  are  caused  by  a  general  reduction  in  the  spatial 
coherence  of  the  original  beam  (Figure  1). 


Coherent 

Phase 

Front 


Turbulence 


Resulting 
Wave  Front 


Figure  1.   Effect  of  Atmospheric  Turbulence  on  Plane  Wave 


B.      STRUCTURE  FUNCTIONS  AND  PARAMETERS 

The  mean  square  difference  between  two  arbitrary  measurements  X  such 
as  temperature  is  defined  as  a  structure  function 

Dx(T1,T2)=<(X(rl)  -x{T2)f>,  (2) 

where  Dx  depends  on  the  separation  distance  [Ref.  8:p.  19].  If  we  assume 
homogeneous,  isotropic  turbulence,  then  Kolmogorov  and  Tatarski  [Ref.  8:p.  46] 
have  shown  that 


dx=cx2i2/\  (3) 


where 


r=\r1-r2\=r12.  (4) 


This  r273  dependence  depicts  Kolmogorov  turbulence  represented  by  the  structure 
parameter  Cx2. 

The  index  of  refraction  structure  parameter,  Cn2  describes  the  fluctuations 
in  the  index  of  refraction  over  a  given  distance.  It  is  the  mean-square  statistical 
average  of  the  difference  in  the  index  of  refraction,  n,  between  two  points 
separated  by  a  distance  r12  [Ref.  9: p.  397] 


L\2 


Similarly,  the  temperature  structure  parameter,  CT2  is 


,      <(T,-T2)2> 
Ct= -T7-, i  (O 


r     2/3 

x12 


where  T,  and  T2  are  the  ambient  atmospheric  temperatures  measured  along  a 
horizontal  path  of  length  r.  According  to  Kolmogorov,  this  relation  is  only  valid 
in  the  inertial  sub  range  described  by  l0<<r<<L0,  where  l0  and  Lg  depend  on  the 
kinetic  energy  and  boundary  conditions.  The  inner  scale,  l0,  where  energy 
dissipation  by  viscosity  and  diffusion  takes  place,  is  typically  on  the  order  of 
millimeters  near  the  ground  and  perhaps  centimeters  in  the  troposphere.  The 
outer  scale,  Lq,  the  largest  scale  of  isotropy,  is  usually  on  the  order  of  meters  or 
10s  of  meters  [Ref.  1:p.  6-9]  and  depends  on  the  boundary  conditions. 

The  index  of  refraction,  n,  is  a  function  of  atmospheric  pressure,  humidity, 
temperature,  and  optical  wavelength.  We  may  neglect  the  effects  of  humidity  in 
all  but  the  most  humid  atmospheric  environments  and  can  disregard  the 
wavelength  dependence  at  optical  wavelengths.  For  optical/IR  wavelengths,  a 
simplified  form  of  this  critical  parameter  is 


,-6P 


72-1  =7  9x1  CT6  —  ,  (7) 

T 


where  P  is  atmospheric  pressure  (millibars)  and  T  is  the  atmospheric  temperature 
(Kelvin)  [Ref.  10:p.  232].  Taking  the  partial  derivative  of  Eq.  (7)  and  assuming 
isobaric  turbulence,  Cn2  depends  on  CT2  by 

Cn2=(79xl0"6-^)2CT2.  (8) 

The  effects  that  turbulence  produces  on  a  laser  beam  are  seen  by 
measuring  the  amplitude  or  phase  fluctuations  of  the  final  image  at  the  target. 
A  common  measure  of  image  quality  is  the  Strehl  ratio  which  is  the  relative 
intensity  of  the  point  spread  function  measured  at  the  mean  centroid  location 
relative  to  an  ideal,  diffraction  limited  optical  system.  The  Strehl  ratio  depends 
on  the  spatial  coherence  of  the  laser  radiation  after  passing  through  the 
atmosphere  which  in  turn  depends  on  the  atmospheric  modulation  transfer 
function.  The  atmospheric  modulation  transfer  function  (MTF)  characterized  by 
the  transverse  coherent  path  length,  r0,  is  a  measure  of  the  electric  field 
transverse  autocorrelation  length,  the  distance  over  which  the  electromagnetic 
field  remains  coherent  [Ref.  1:p.  6-29].    Fried  showed  that,  for  plane  waves, 

r0=2.1[1.46^2fcn2(z)dz]-3/5,  (9) 

o 

where  k  is  the  optical  wave  number  (2rr/A)  and  L  is  the  optical  path  length  [Ref. 
9: p.  397].     Typical  values  of  r0  range  from  centimeters  in  regions  of  high 
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turbulence,  to  tens  of  centimeters  at  good  astronomical  sites,  to  meters  and 
perhaps  larger  in  low  turbulence  regions. 

C.      CURRENT  APPLICATIONS 

Some  alternatives  to  the  TMD  scenario  involve  near  horizontal  atmospheric 
propagation  at  altitudes  of  10-20  km.  One  initiative,  the  Airborne  Laser  Program 
(ABL),  involves  mounting  a  laser  in  an  aircraft  flown  above  the  tropopause  into 
the  stratosphere  to  avoid  clouds,  sub-visual  cirrus  (ice  crystals),  and  tropospheric 
turbulence  [Ref.  11].  The  upper  altitude  limit  arises  from  aircraft  altitude 
constraints  although  other  platforms  such  as  low  orbit  satellites  are  also  being 
considered  (space-based  alternatives  are  subject  to  more  restrictive  political 
constraints).  The  lower  altitude  limit  essentially  is  the  bottom  of  the  stratosphere. 
The  reason  for  this  is  that  beam  propagation  in  the  troposphere  encounters 
higher  turbulence,  particularly  near  the  jet  stream,  and  has  higher  absorption  and 
scattering  from  atmospheric  aerosols,  water  vapor,  ice  crystals,  and  clouds. 

With  adequate  phase  information,  adaptive  optics  can  pre-distort  the 
wavefront  of  the  beam  to  correct  for  turbulence  that  is  along  the  propagation 
path.  This  is  necessary  to  achieve  a  30  cm  -  1.0  m  diameter  spot  size  on  a 
target  at  a  range  of  400  km  (0.8  /yrad).  The  Airborne  Laser  Program  needs  an 
accurate  assessment  of  the  atmospheric  turbulence  along  stratospheric  optical 
paths  to  complete  a  systems  analysis.  This  was  the  driving  force  for  this 
investigation. 


D.       PREVIOUS  RESULTS 

GP  and  NPS  have  conducted  numerous  experiments  to  measure  the 
temperature  structure  parameter  using  balloon-borne  thermosondes  and  single- 
point  temperature  probes.  Thermosondes  measure  differential  temperature 
fluctuations  across  a  horizontal  path,  commonly  one  meter  [Ref.  12:p.  14].  The 
mean  square  fluctuation  in  temperature  between  two  3.8  /ym  platinum-coated, 
tungsten  probes  separated  by  one  meter  gives  a  measurement  of  CT2  which,  in 
turn,  gives  Cn2  using  Eq.  (8)  [Ref.  13:p.  2].  The  NPS  single-point  temperature 
probes  measure  temperature  using  100/ym  thermistors  along  both  the  ascending 
and  descending  trajectories  of  the  balloon. 

Previous  experiments  and  analyses  have  found  a  clear  increase  in  the  index 
of  refraction  structure  parameter  Cn2  above  -8  km  during  daylight  hours.  In 
June  1990,  GP  conducted  a  series  of  thermosonde  launches,  in  conjunction  with 
NPS  optical  measurements.  Figures  2  and  3  show  plots  of  Cn2  from  both  day 
and  night  launches.  The  raw  data  minima  in  Figure  2,  a  night  flight,  approach 
the  sonde  noise  level  (solid  curve)  for  much  of  the  flight.  The  raw  data  in  the 
daytime  flight  (Figure  3),  almost  never  approach  the  sonde  noise  level  [Ref.  3:p. 
55].    This  anomaly  has  been  characteristic  of  all  day  flights  and  is  suspicious. 
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Figure  2.   Log  (Raw  Cn2)  Profile  Collected  by  the  Geophysics 
Directorate  on  Maui,  2  June  1990  (Night)  [Ref.  3:p.  55]. 
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Figure  3.   Log  (Raw  Cn2)  Profile  Collected  by  the  Geophysics 
Directorate  on  Maui,  3  June  1990  (Day)  [Ref.  3:p.55]. 
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GP  and  Sparta,  Inc.,  published  a  scientific  report  in  August  1990  addressing 
temperature  perturbations  from  solar  heating  of  thin  layers  of  aerosol  and  ozone. 
They  found  that  there  was  some  evidence  to  support  the  postulate  of  temperature 
perturbations  due  to  solar  heating  of  ozone,  although  this  could  not  be  proven 
due  to  the  low  vertical  resolution  nature  of  03.  GP  concluded  that  it  was  not 
probable  that  variations  in  CT2  were  due  to  perturbations  in  background 
stratospheric  aerosols.  [Ref.  13:p.  14]  GP  also  investigated  the  possibility  that 
solar  radiation  heating  of  the  probe  support  wires  caused  instrumental  effects  but 
concluded  that  this  was  not  responsible  for  the  prevalent  diurnal  variation  [Ref. 
2:p.  55]. 

Previous  work  at  NPS  suggested  that  solar  heating  and  probe  wire 
vibrations,  caused  by  the  wind  from  the  balloon's  ascent  coupled  with  the  probe 
package's  oscillating  motion,  produced  an  AC  signal  that  seriously  degraded  the 
daytime  data  [Ref.  3].  An  actual  increase  in  daytime  turbulence  was  extremely 
unlikely  since  the  day-night  change  in  turbulence  occurred  within  minutes  of 
sunrise  or  sunset.  It  is  unlikely  that  this  rapid  onset  and  decay  is  a  real 
phenomena  because  large-scale  molecular  diffusion  of  energy  throughout  the 
atmosphere  should  not  be  possible  on  such  small  time  scales.  Optical 
measurements  collected  in  July  of  1990  showed  no  indication  of  a  change  in 
isoplanatic  angles  between  day  and  night  whereas  thermosonde  results  collected 
during  the  same  period  showed  persistent  day-night  differences.  Optical 
isoplanatic  angles  are  sensitive  to  Cn2  weighted  by  altitude  z5/3  [Ref.  3:p.  9]. 
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Consequently,  turbulence  in  the  5-15  km  region  has  the  dominant  contribution. 
Weitekamp  showed  that  simulated  solar  heating  of  a  GP  3.8  /ym  tungsten  probe 
in  the  laboratory  produced  a  0.1  to  0.2  C  temperature  increase.  This  was 
sufficient  to  account  for  the  observed  increase  in  the  day  CT2  values  when 
coupled  with  probe  vibration  and  sonde  rotary  and  pendulum  motion. 
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III.    DIURNAL  VARIATION 

A.      DIURNAL  VARIATION  EXPERIMENT  2 

1.      Purpose 

The  Geophysics  Directorate  group  conducted  a  series  of  thermosonde 
launches  from  4  to  9  December  1990  at  Holloman  Air  Force  Base,  New  Mexico, 
to  study  the  diurnal  variation  in  turbulence  measurements  of  Cn2.  The  purpose 
was  to  validate  or  disprove  the  existence  of  measurement  artifacts  that  were 
responsible  for  the  diurnal  variations  observed  by  both  NPS  and  GP  in  previous 
experiments.  [Ref.  14]   These  artifacts  include: 


1.  Balloon  thermal  wake  in  combination  with  payload  motion  (pendulum 
and  rotation). 

2.  Differential  solar  heating  of  the  probes  due  to  misalignment  in 
combination  with  payload  motion  (pendulum  and  rotation). 

3.  Differential  solar  heating  of  the  probes  due  to  different  wire  geometry  or 
emissivity/reflectivity  in  combination  with  payload  motion  (pendulum  and 
rotation). 

4.  Differential  solar  heating  of  the  probes  due  to  probe  wire  oscillation, 
possibly  exacerbated  by  payload  motion  (pendulum  and  rotation). 


2.       Experimental  Procedure 

The  GP  thermosonde  system  measured  the  temperature  difference 
between  two  3.8  //m  platinum-coated  tungsten  wire  probes  exposed  to  the 
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ambient  atmosphere.  The  probe  wires  were  flexible  s-shaped  or  u-shaped  wires 
approximately  0.45  cm  long  connected  to  two  probe  supports.  The  probes 
formed  elements  of  an  AC  Wheatstone  bridge  which  measured  the  difference  in 
resistance  between  them.  High  and  low  pass  filters  restricted  the  system 
bandwidth  to  the  range  .2  Hz  to  1000  Hz  at  the  -3dB  level.  The  thermosonde 
passed  the  temperature  difference  T2  -  T,  through  a  root  mean  square  (RMS) 
module  with  a  4-second  time  constant  to  ground-based  monitoring  equipment 
every  four  seconds  which  logged  the  data  and  calculated  Cn2  [Ref.  12:p.  15]. 
The  sondes  ascended  at  an  average  rate  of  approximately  5  meters  per  second. 
An  ascent  to  30  km  required  approximately  95  minutes.  Historically,  a  GP 
thermosonde  launch  consisted  of  a  thermosonde  package  with  a  single  pair  of 
probes  with  a  1  m  spacing  connected  to  a  VIZ  MICROSONDE  1  telemetry 
package.  The  December  1990  launches  were  unique  since  they  used  two 
separate  thermosonde  packages  strapped  together  with  a  single  VIZ  telemetry 
package.  The  two  pairs  of  probes  were  arranged  in  1  m  differential  pairs  on  a 
single  50  x  50  mm  styrofoam  boom,  so  that  each  sensor  element  was 
approximately  8  cm  from  the  corresponding  member  of  the  other  pair.  There 
were  two  probe  configurations  for  the  experiment.  On  most  flights,  both  pairs  of 
probes  were  parallel.  One  flight  had  one  pair  of  probes  parallel  and  the  second 
set  perpendicular  to  each  other. 
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3.      Experimental  Conditions 

In  all,  five  successful  balloon  launches  occurred  during  December 
1990.  Four  were  day  launches  and  one  was  a  night  launch.  Local  sunset 
occurred  around  1653  hours.  Difficulties  in  controlling  the  heavier  dropsonde 
packages  precluded  any  data  taking  during  descent.  Table  I  contains  a 
summary  of  the  individual  sonde  launches. 

a.  4  December  1990 

Time  of  launch  of  DVE9465  was  1450  hours  local  and  both  sets 
of  probes  were  parallel.   The  last  good  data  were  obtained  at  29.03  km. 

b.  5  December  1990 

The  time  of  launch  of  DVE9460  was  1407  hours  local  on  5 
December  1990  and  both  sets  of  probes  were  parallel.  Good  data  were  obtained 
up  to  26.48  km. 

c.  6  December  1990 

The  time  of  launch  of  DVE9454  was  1330  hours  local  on  6 
December  1990.  Both  sets  of  probes  were  oriented  in  a  parallel  fashion.  The 
maximum  altitude  providing  useful  data  was  20.5  km. 

d.  7  December  1990 

The  time  of  launch  of  DVE9448  was  1357  hours  local  on  7 
December  1990.    One  set  of  probes  was  in  a  parallel  configuration  while  the 
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second  set  was  perpendicular.    This  was  the  only  launch  in  which  this  was  the 
case.    The  maximum  altitude  that  provided  useful  data  was  28.86  km. 

e.    8  December  1990 

DVE9462  was  the  only  night  launch.  Time  of  launch  was  1837 
hours  local.  Both  sets  of  probes  were  parallel.  Good  data  were  taken  up  to 
30.02  km. 


TABLE  I.  Diurnal  Variation  Experiment  2  Launch  Summary,  Conducted  by  the 
Geophysics  Directorate  at  Holloman,  Air  Force  Base,  NM  from  4-8  December 
1990. 


DATE 

4DEC  90 

5DEC  90 

6DEC  90 

7DEC  90 

8DEC90 

FLIGHT 

DVE9465 

DVE9460 

DVE9454 

DVE9448 

DVE9462 

TIME 
(Local) 

1450 

1407 

1330 

1357 

1837 

PROBES 

par/par 

par/par 

par/par 

par/perp 

par/par 

MAX  ALT 
(km) 

29.03 

26.48 

20.5 

28.86 

30.02 
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B.       RESULTS  AND  ANALYSIS  OF  DIURNAL  VARIATION  EXPERIMENT  2 

1.  Data  Reduction 

The  telemetry  package  continuously  transmitted  data,  at  four  second 
intervals,  to  the  ground-based  monitoring  equipment,  resulting  in  data  sets  of 
approximately  1450  lines  for  a  90  minute  flight.  The  data  were  a  succession  of 
independent  4  second,  RMS  temperature  differences  in  the  form 

[<(r2-r1)2>]1/2  (io) 

Eq.  (6),  the  temperature  structure  function  equation,  along  with  the  temperature 
difference  data  and  the  r=1  m  spacing  of  the  probes  gave  a  CT2  value  for  each 
altitude  at  which  measurements  were  made,  approximately  every  20  m. 
Combining  these  data  with  measured  pressure  and  temperature  yielded  Cn2 
values  at  each  altitude. 

2.  Statistical  Methodology 

To  determine  if  any  measurement  artifacts  contaminated  the  data,  we 
performed  three  statistical  analyses  of  the  data  taken  from  the  two  pairs  of 
probes.  The  first  involved  examining  plots  of  the  individual  data  points  vs. 
altitude.  This  can  reveal  major  inconsistencies.  The  second  involved  an 
examination  of  x-y  cross-correlation  plots  between  data  from  pairs  of  probes. 
The  third  involved  computation  of  the  correlation  coefficient  between  probe  pairs. 
The  temperature  differences  measured  at  each  successive  altitude  should  have 
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a  high  correlation,  since  both  sets  of  probes  were  essentially  passing  through  the 
same  air  mass.  If  there  were  independent  measurement  artifacts,  however; 
inconsistencies  in  the  x-y  cross-correlation  plots  and  the  correlation  between  data 
could  possibly  provide  some  insight  as  to  the  source  of  the  error. 

This  thesis  used  the  Microsoft  EXCEL  spreadsheet  to  compute  the 
statistical  correlation  between  both  pairs  of  RMS  temperature  fluctuations.  The 
correlation  coefficient,  r,  was  computed  according  to  the  standard  statistical 
expression 


I{CTX,CT2)  =    <3cT1-CT2   ,  (n) 


\l°cri° 


CT2 


where 


%'^-iE<CT1-CT1)(CPl-Cra)f  (12) 


and  CT1  was  the  RMS  temperature  difference  for  one  pair  of  probes  and  CT2 
was  the  difference  for  the  second  pair.  Values  of  r~  1  indicate  a  high  degree  of 
correlation;  values  of  r~0  indicate  little  or  no  correlation. 

For  analyses,  the  30  km  data  sets  were  divided  into  5  km  segments. 
We  also  computed  a  single  r  coefficient  for  each  flight  that  contained  all  the  data 
points.  Applying  Eq.  (12)  to  each  vertical  5  km  region  (or  to  the  complete  flight 
in  the  case  of  the  total  data  set)  gave  the  standard  deviation  needed  in  Eq.  (11), 
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which    in    turn,    yielded    the    correlation    coefficient    between    each    pair    of 
measurements. 

The  measured  noise  level  of  the  CT2  instrument  packages  was  .002K. 
Consequently,  (AT)^  is  undefined  for  values  less  than  .002K  [Ref.  12:p.  14].  All 
data  taken  that  produced  a  reading  of  <.002K  were  discarded  as  being 
statistically  invalid  because  it  would  not  be  possible  to  discriminate  between  the 
signal  and  the  noise. 

3.       Results 

a.    Diurnal  Variation 

Plots  of  Log  (Raw  Cn2)  vs.  altitude  show  the  typical  order  of 
magnitude  increase  in  optical  turbulence  in  the  day  compared  to  night  in  the  10- 
20  km  region.  Figures  4  and  5  of  the  7  December  flight  represent  a  daytime 
flight  and  are  also  the  only  plots  obtained  using  the  parallel/perpendicular 
configuration  (Appendix  A  contains  the  other  daytime  plots).  Figures  6  and  7 
show  the  data  taken  by  the  two  sets  of  probes  on  the  8  December  night  flight. 
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Figure  4.  Log  (Raw  Cn2)  vs.  Altitude  Profile  (Daytime)  , 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  30  km, 
Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set  1) . 
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Figure  5.  Log  (Raw  Cn2)  vs.  Altitude  Profile  (Daytime), 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  3  0  km, 
Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set  2) . 
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Figure  6.  Log  (Raw  Cn2)  vs.  Altitude  Profile  (Night), 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  30  km, 
Holloman  Air  Force  Base,  NM,  8  December  1990  (Probe  Set  1) . 
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Figure  7.  Log  (Raw  Cn2)  vs.  Altitude  Profile  (Night), 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  3  0  km, 
Holloman  Air  Force  Base,  NM,  8  December  1990  (Probe  Set  2) . 
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b.    7  December  1990 

On  7  December  1990,  a  day  flight  with  both  parallel  and 
perpendicular  probe  configurations,  the  two  different  sets  of  probes  behaved 
differently.  The  set  of  parallel  probes  exhibited  a  trend  common  to  the  other  day 
flights.  At  low  altitudes,  up  to  approximately  8-12  km,  some  of  the  data  points 
closely  approached  the  baseline  temperature-pressure  curve.  As  the  probe 
continued  to  ascend,  the  log  (Raw  Cn2)  values  increased  and  moved  away  from 
the  system  noise  level  curve  by  an  order  of  magnitude  (Figure  4).  The 
perpendicular  set  of  probes  exhibited  enhanced  differences.  These  log  (Raw  Cn2) 
values  consistently  were  an  order  of  magnitude  higher  than  the  noise  level 
throughout  the  entire  flight  (Figure  5).  This  indicates  that  there  was  a  constant 
source  of  error,  or  contamination,  throughout  the  entire  flight. 

Figures  8  thru  1 1  show  the  x-y  correlation  between  the  two  values 
(CT1  and  CT2)  measured  at  each  altitude.  These  figures  show  the  effect  of  noise 
or  some  constant  signal  with  a  value  of  approximately  0.006  K  that  affects  the 
data  in  channel  2.  Examination  of  1 0  km  increments  of  the  atmosphere  reveals 
that  this  signal  first  became  apparent  between  10-20  km  (Figure  10)  and  became 
very  pronounced  at  20-30  km  (Figure  11).  The  dense  spike  in  the  data  at 
CT1  «  0.006  K  illustrates  this  effect. 
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CT1  vs.  CT2  (DAY) 

Par/Perp  Configuration  (7  DEC  1990) 
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Figure  8.  CT1  vs.  CT2  (Day),  Correlation  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 
(Parallel/Perpendicular  Configuration)  Collected  Between  1.2 
and  30  km,  Holloman  Air  Force  Base,  NM,  7  December  1990. 
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CT1  vs.  CT2  (DAY) 

Par/Perp  Configuration  (7  DEC  1990) 
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Figure   9.      CT1   vs.    CT2    (Day),    Correlation   of   RMS   Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 

(Parallel/Perpendicular   Configuration)    Collected   Between    1.2 
and    10   km,    Holloman  Air   Force   Base,    NM,    7    December    1990. 
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CT1  vs.  CT2  (DAY) 

Par/Perp  Configuration  (7  DEC  1990) 
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Figure   10.      CT1  vs.    CT2    (Day),    Correlation  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 

(Parallel/Perpendicular    Configuration)     Collected    Between    10 
and   20   km,    Holloman  Air   Force   Base,    NM,    7   December    1990. 
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CT1  vs.  CT2  (DAY) 

Par/Perp  Configuraiion  (7  DEC  1990) 
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Figure   11.      CT1  vs.    CT2    (Day),    Correlation  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 

(Parallel/Perpendicular    Configuration)     Collected    Between    20 
and    30   km,    Holloman  Air   Force   Base,    NM,    7    December    1990. 
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The  correlation  coefficient  for  the  entire  flight  was  50%,  the  lowest 
of  the  five  successful  flights.  Table  II  shows  that  r  was  high,  96%,  in  the  first  5 
km  interval  and  was  below  50%  for  the  rest  of  the  flight. 


Table  II.  CT1  vs.  CT2  Correlation  for  7  December  1990  Daytime  Launch  with 
Parallel/Perpendicular  Probe  Configuration  at  Holloman  Air  Force  Base,  NM. 


ALTITUDE  (km) 

CORRELATION,    r 

Data  points 

1.2-5 

.96 

155 

5-10 

.42 

220 

10-15 

.32 

252 

15-20 

.27 

174 

20-25 

.34 

240 

25-30 

.40 

176 

TOTAL 

.50 

1217 

c.    8  December  1990 

The  plots  of  log  (Raw  Cn2)  values  for  the  8  December  1990  flight 
followed  the  characteristic  nighttime  behavior.  Some  of  the  data  consistently 
approached  the  noise  level  up  to  approximately  15-20  km,  at  which  point  they 
showed  a  slight  shift  to  the  higher  values.   Figures  6  and  7  show  these  results. 
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Figures  12  thru  15  plot  the  CT1  vs.  CT2  correlation  vs.  altitude. 
These  figures  show  the  effect  of  a  small  but  relatively  constant  offset.  The  CT2 
values  are  consistently  greater  than  the  CT1  values.  This  could  not  occur  without 
a  signal,  instrument,  or  calibration  error  since  the  probes  passed  through 
essentially  the  same  air  mass  up  to  30  km  and  took  893  independent 
measurements.  Table  III  shows  that  the  correlation  coefficient  was  high  in  all 
intervals  with  the  exception  of  the  5-10  km  region  which  had  an  insufficient 
number  of  data  points  above  the  noise. 


TABLE  III.  CT1  vs.  CT2  Correlation  for  Geophysics  Directorate  8  December  1990 
Night  Launch  with  Dual-Parallel  Probe  Configuration  at  Holloman  Air  Force  Base, 
NM. 


ALTITUDE  (km) 

CORRELATION,  r 

Data  Points 

1.2-5 

.98 

115 

5-10 

.58 

29  * 

10-15 

.90 

161 

15-20 

.92 

240 

20-25 

.86 

199 

25-30 

.86 

148 

TOTAL 

.91 

893 

statistically  insufficient 
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CT1  vs.  CT2  (NIGHT) 

Dual-Parallel  Config.  (8  DEC  1990) 
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Figure  12.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperature  Fluctuations 
Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration)  Collected  Between 
1.2  and  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1990. 
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CT1  vs.  CT2  (NIGHT) 

Dual-Parallel  Config.(8  DEC  1990 
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Figure  13.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperature  Fluctuations 
Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration)  Collected  Between 
1.2  and  10  km,  Holloman  Air  Force  Base,  NM,  8  December  1990. 
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CT1  vs.  CT2  (NIGHT) 

Dual-Parallel  Config.  (8  DEC  1990) 
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Figure  14.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperature  Fluctuations 
Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration)  Collected  Between 
10  and  20  km,  Holloman  Air  Force  Base,  NM,  8  December  1990. 
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CT1  vs.  CT2  (NIGHT) 

Dual-Parallel  Config.  (8  DEC  1990 
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Figure  15.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperature  Fluctuations 
Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration)  Collected  Between 
20  and  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1990. 
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The  plots  of  the  4,  5,  and  6  December  flights  followed  the  typical 
daytime  trend  and  are  in  Appendix  A.  The  data  collection  during  the  descent 
phase  of  the  experiment  was  not  successful  due  to  difficulties  in  using  the  heavier 
than  standard  package. 

4.      Analysis 

The  high  degree  of  correlation  of  all  flights  in  the  1-5  km  region  is  due 
in  part  to  the  turbulent  boundary  layer  which  produces  high  Cn2  in  the  daytime. 
Following  Weitekamp,  solar  induced  temperature  differences  due  to  insolation  are 
smallest  at  ground  level  and  are  the  largest  at  the  tropopause  when  the  probe 
is  at  a  temperature  of  approximately  200  K.  During  the  day,  solar  energy  raises 
the  probe  surface  temperature  above  the  ambient  air  temperature.  This 
difference  increases  as  the  probe's  temperature  declines.  This  increase  mimics 
turbulent  fluctuations  when  it  varies  rapidly  with  time.  At  night,  however,  the 
surface  temperature  drops  below  that  of  the  air  as  it  cools  by  radiation.  The 
decrease  in  correlation  in  the  5-10  km  region  appears  to  be  the  result  of  the 
decline  in  the  boundary  layer  contribution  coupled  with  the  increase  in  the  solar 
induced  probe-atmosphere  temperature  differences. 

The  only  flight  with  the  parallel/perpendicular  configuration  of  the  two 
sets  of  probes,  launched  on  7  December,  had  by  far  the  lowest  correlation  (total 
average  r  of  50  %).  This  indicates  that  the  probe  orientation  with  respect  to  the 
sun  influenced  the  differences  in  the  measured  temperatures  [Ref.  3].  The 
motion  of  the  ascending  sonde  package  continuously  changed  the  orientation, 
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and  hence  the  projected  surface  area,  of  the  probes  with  respect  to  the  incident 
solar  radiation  resulting  in  uneven  thermodynamic  heating  and  cooling  of  the 
probes  [Ref.  3:p.  33].  The  results  of  the  7  December  1990  flight  further  support 
his  findings. 

The  most  significant  results  occurred  in  the  region  of  most  pressing 
interest,  namely  10  to  20  km.  The  nighttime  data  correlation  was  approximately 
91  %,  markedly  higher  than  any  of  the  daytime  data  (Figures  10  and  14).  This 
indicated  that  solar  contamination  of  the  daytime  data  is  a  real  effect. 

The  4  December  1990  flight  had  very  low  correlation  in  the  15  km  to 
25  km  region;  much  lower  than  even  the  7  December  flight  (see  Table  IV).  This 
could  have  been  due  to  instrument  error,  calibration,  or  perhaps  the  probes 
actually  measured  an  increase  in  turbulence.  Some  sort  of  instrument 
malfunction  was  most  likely  the  case.  The  fact  that  the  correlation  of  this  flight 
jumped  to  77%  from  25  km  to  30  km,  compared  to  16%  in  the  20  to  25  km 
altitude  region  indicates  a  systematic  error.  This  marked  increase  was  the  most 
dramatic  data  discontinuity  between  any  altitude  regions  in  any  of  the  five  flights 
in  Diurnal  Variation  Experiment  2  and  should  be  considered  skeptically  prior  to 
making  any  inferences  from  it. 

Figure  1 6  is  a  plot  of  the  results  of  the  correlation  analysis  of  each 
flight.   Table  IV  summarizes  the  data  for  the  4,  5,  and  6  December  flights. 
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CORRELATION  OF  CT1  AND  CT2 

r  vs.  5  km  Vertical  Regions 

7  DEC  day  flight  was  only  flight  with  parallel/perpendicular  configuration. 
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Figure  16.  Correlation  of  CT1  and  CT2  vs.  Altitude  for  all  Diurnal  Variation 
Experiment  Flights  Conducted  by  the  Geophysics  Directorate  at  Holloman,  Air 
Force  Base,  NM,  from  4-8  December  1990. 
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TABLE  IV.    CT1  and  CT2  Correlation  vs.  Altitude  for  Probe  Pairs  for  Remaining 
Flights  in  Diurnal  Variation  Experiment  2. 


DATE 

ALT  (km) 

CORRELATION  (r) 

#  of  POINTS 

4  DEC  90 

1.2-5 

.98 

148 

5-10 

.79 

52 

10-15 

.58 

262 

15-20 

.19 

241 

20-25 

.16 

229 

25-30 

.77 

165 

TOTAL 

.63 

1285 

5  DEC  90 

1.2-5 

.71 

11 

5-10 

.79 

52 

10-15 

.39 

155 

15-20 

.73 

154 

20-25 

.71 

107 

25-30 

.60 

9 

TOTAL 

.82 

488 

6  DEC  90 

1.2-5 

.97 

232 

5-10 

.62 

217 

10-15 

.77 

60 

15-20 

.54 

41 

TOTAL 

.78 

550 
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C.      DVE  SUMMARY 

The  high  correlation  obtained  for  the  temperature  differences  (CT1  and  CT2) 
for  the  sets  of  probes  in  the  night-time  flight  on  8  December  1990  and  the  low 
correlation  obtained  during  day-time  flights  particularly  near  the  tropopause 
support  the  hypothesis  that  the  apparent  order  of  magnitude  increase  in  optical 
turbulence  is  an  instrument  effect  attributable  to  solar  heating  of  the  probes.  This 
is  consistent  with  previous  work  that  shows  that  the  solar  induced,  probe-air 
temperature  differences  are  largest  at  low  temperatures  [Ref.  3:p.  19].  The  fact 
that  the  only  flight  that  had  the  parallel/perpendicular  configuration  of  the  two 
sets  of  probes,  the  7  December  launch,  had  the  largest  overall  variations 
(smallest  correlation  coefficient,  r),  indicates  that  the  probe  orientation  with 
respect  to  the  sun  contributed  to  the  differences  in  the  measured  temperatures. 

The  motion  of  the  package  (rotation  or  pendulum-like  oscillations),  could 
produce  a  low  frequency  signal.  The  parallel/perpendicular  configuration  would 
correlate  more  poorly  than  the  double  parallel  configuration  because  the  probes 
would  move  through  the  air  masses  with  different  orientations  and  be  exposed 
to  different  conditions.  Differential  solar  heating  of  the  probes  from  probe  wire 
oscillation,  possibly  exacerbated  by  payload  motion  (pendulum  and  rotation) 
could  produce  a  higher  frequency  signal  which  could  mask  the  data  [Ref.  3: pp. 
16,46].  The  flexible  probe  wires,  which  oscillate  from  the  cross  flow  of  air,  would 
continuously  present  differing  cross-sectional  areas  to  the  sun  and  induce  a 
signal.  If  probe  vibrations  were  the  cause  of  these  differing  results,  then  both  the 
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parallel/perpendicular  configuration  and  the  dual  parallel  configuration  would 
have  exhibited  uneven  heating  behavior.  The  daytime  data  currently  available  are 
only  reliable,  meaning  consistent,  with  nighttime  data  up  to  approximately  8-10 
km,  which  is  an  altitude  that  is  clearly  beneath  the  tropopause.  Daytime 
microthermal  probe  data  beyond  this  are  highly  questionable.  Although  the 
number  of  flights  in  this  series  of  launches  was  limited,  these  results  are 
consistent  with  those  of  many  previous  day  launches  and  support  the  hypothesis 
that  solar  heating  of  the  probes  contaminates  the  daytime  data.  Although  the 
exact  mechanism  of  this  contamination  is  not  readily  apparent,  the  evidence 
confirming  its  existence  is  convincing. 

If  balloon  wake  (thermal  shedding)  was  responsible  for  the  diurnal  variation 
in  the  CT2  measurements,  the  correlation  of  the  nighttime  flights  should  be  poorer 
than  the  day.  The  balloon,  which  is  cooler  than  its  surroundings  at  night,  creates 
a  cooled  wake  beneath  it  as  it  rises.  The  wake  effect  should  be  worse  at  night 
since  the  temperature  difference  between  the  balloon  and  the  surrounding  air  is 
larger  than  in  the  day.  Weitekamp  showed  that  this  effect  was  small  at  the 
distances  that  the  thermosonde  packages  were  suspended  beneath  the  balloon 
[Ref.  3:p.  25]. 

D.      SONDE  STABILIZATION  FOR  SUBSEQUENT  EXPERIMENTS 

To  eliminate  solar  effects  apparent  in  the  daytime  launches  due  to  rotation 
and  oscillation  of  the  sonde,  it  is  necessary  to  stabilize  the  sonde  as  much  as 
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possible.  The  most  reliable  and  consistent  data  will  occur  when  the  probes 
maintain  a  consistent  alignment  with  respect  to  the  sun.  Possible  solutions  to 
this  may  include  the  use  of  a  gyro  mechanism  to  stabilize  the  platform,  vertical 
baffles  or  fin  stabilizers,  solar  cells  to  determine  relative  solar  angle,  and 
weights/anchors  to  hold  the  apparatus  steady.  To  reduce  the  effects  of  probe 
vibration,  the  tungsten  wire  probes  must  be  attached  straight  across  the  probe 
supports  rather  than  with  a  flexible  configuration.  This  will  involve  a  delicate 
compromise  to  avoid  thermal  stress  and  strain  on  the  wire  and  its  contacts  and 
yet  suppress  vibration. 
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IV.    PROBE  RESPONSE 

A.      PROBE     RESPONSE    TIME    AS    A    FUNCTION     OF    DECREASING 
TEMPERATURE  AND  AIR  DENSITY 

1.      Thermocouple  Measurement  Limitations  Due  to  Probe  Response 
Time 

A  second  systematic  issue  that  impacts  the  performance  of  the  thermal 
probe  turbulence  measurements  involves  estimating  the  probe  response  time  vs. 
altitude.  When  a  thermocouple  measures  the  temperature  of  its  atmospheric 
surroundings,  the  value  that  it  provides  is  its  own  temperature,  not  that  of  the 
surrounding  air.  This  occurs  because  of  the  finite  heat  transfer  mechanisms 
occurring  at  the  surface-air  interface.  The  thermal  mass  of  the  probe  and  the 
convective  air-surface  interface  limit  the  rate  of  heat  transfer  producing  a  delayed 
probe  response.  This  limits  the  high  frequency  components  of  the  signal  used 
to  compute  structure  functions.  [Ref.  15:p.  15]  This  becomes  increasingly 
important  for  structure  function  measurements  with  small  separation  distances. 
A  10  cm  separation  requires  a  probe  response  time  of  1-2  ms.  This  problem 
increases  at  higher  altitudes.  The  low  density  environment  in  the  10  to  20  km  (10 
mb)  region  offers  less  effective  heat  transport  (primarily  forced  convection)  by 
which  the  probes  adjust  to  the  temperature  of  their  surroundings,  consequently 
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measurements  underestimate  the  corresponding  structure  functions  as  altitudes 
increase. 

2.       Fluid  Characteristics  and  Heat  Transfer 

a.    Forced  Convection 

Forced  convection  is  the  dominant  (albeit,  still  a  weak)  heat 
transfer  mechanism  in  regions  where  the  mean  free  path  of  the  surrounding 
molecules  is  comparable  to  the  dimensions  of  the  heat  transfer  surface  (the 
cylindrical  wire  probe).  The  physical  mechanism  of  forced  convection  which 
removes  heat  (energy)  from  the  probes  is  the  collisional  momentum  transfer  to 
surrounding  air  molecules  when  these  molecules  impact  the  probe  [Ref.  16:p. 
631].  The  mean  free  path  of  atmospheric  molecules  increases  from  0.06  //m  at 
sea  level  to  approximately  4.4  /vm  at  30  km  [Ref.  17:p.  14-13].  This  is 
comparable  to  the  dimensions  of  a  12.4  /ym  diameter  cylindrical  wire  probe. 
These  are  the  dimensions  of  the  ANSI  Type  E  chromel-constantan  thermocouples 
that  the  NPS  group  used  in  the  summer  and  fall  1992  experiments;  the 
thermocouples  were  manufactured  by  Omega  Engineering,  Inc.  [Ref.  18:p.  A-33]. 
Above  30  km,  as  the  molecular  mean  free  path  length  continues  to  increase, 
forced  convection  is  no  longer  the  dominant  heat  transfer  mechanism  for  probes 
of  this  size. 
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b.    Viscosity 

The  thermodynamic  transport  property  that  characterizes  the  fluid- 
mechanical  behavior  of  the  medium  is  the  viscosity,  jj,  also  known  as  the 
dynamic  viscosity.  It  exists  as  a  point  function  in  a  flowing  fluid,  the  air,  and 
follows  all  of  the  laws  and  state  relations  of  ordinary  equilibrium  thermodynamics 
[Ref.  7:p.  20].  Viscosity  is  a  measure  of  the  diffusion  of  momentum  parallel  to  the 
flow  velocity  and  transverse  to  the  gradient  of  the  flow  velocity  [Ref.  16:p.  219] 
and  has  the  dimensions  of  (kglm>s).  It  relates  the  local  stresses  in  a  moving 
fluid  to  the  rate  of  strain  in  the  fluid  element  by 

T=H^,  (13) 

dy 

where  r  is  the  shear  stress  and  du/dy  is  the  local  fluid  velocity  gradient  [Ref.  7:p. 
26]. 

Prandtl  showed  that  low  viscosity  fluid  flows  such  as  airflow  are 
divided  into  a  thin  viscous  boundary  layer  near  a  solid  surface  which  borders  a 
nearly  inviscid  outer  layer  which  is  still  governed  by  the  Euler  and  Bernoulli 
Equations  [Ref.  19:p.  46]  Under  shearing,  a  fluid  begins  to  move  at  a  strain  rate 
inversely  proportional  to  the  coefficient  of  viscosity.  Viscosity  depends  primarily 
on  temperature  at  low  altitudes  by 

r        T+S 
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where  0  is  a  constant  equal  to  1.458  x  106  kgs"1m"1K'1/2  and  S,  Sutherland's 
constant,  is  equal  to  110.4  K  [Ref.  20:p.  14-5]. 

c.    Kinematic  Viscosity 

The  kinematic  viscosity,  v,  is  the  ratio  of  the  dynamic  viscosity  of 
a  gas  to  its  density  by 


v=A 


(15) 


where  p  is  the  atmospheric  density  in  (kg/m3)  [Ref.  20:p.  14-6].   The  dimensions 
of  v  are  (m2/s).   Figure  17  shows  how  kinematic  viscosity  increases  with  altitude. 
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Figure  17.    Kinematic  Viscosity  as  a  Function  of  Altitude  [Ref.  20:p.  14-6]. 
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of.    No-Slip  Condition 

When  a  fluid  flow  encounters  a  solid  surface,  molecular 
interactions  cause  the  fluid  and  the  surface  to  seek  energy  and  momentum 
equilibrium  with  each  other;  all  liquids  are  essentially  in  equilibrium  with  the 
surface  with  which  they  are  in  contact;  fluids  in  contact  with  a  solid  take  on  the 
velocity  and  temperature  of  the  surface.  These  are  the  no-slip  and  no- 
temperature-jump  conditions  which  are  only  valid  when  the  molecular  mean  free 
path  is  less  than  the  surface  dimensions  of  the  solid.  [Ref.  7:p.  35] 

e.    Reynolds  Number,  Re 

The  parameter  that  describes  the  viscous  behavior  of  a  Newtonian 
fluid  such  as  air  is  the  dimensionless  Reynolds  number  Re 

*-£*.*,  (16) 

\i  V 

where  p  is  the  atmospheric  density,  jj  is  the  dynamic  viscosity,  V  is  the 
characteristic  velocity,  and  L  is  the  characteristic  length  scale  of  the  flow.  The 
Reynolds  number  is  a  ratio  of  the  dynamic  forces  of  mass  flow  (inertia)  to  the 
viscous  forces  within  the  atmosphere.  For  a  12.4  //m  chromel-constantan  probe, 
the  value  of  the  Reynolds  number  varies  from  3.75  at  sea  level  to  0.07  at  30  km 
for  windspeeds  of  4.4  m/s.  Figure  18  is  a  plot  of  the  Reynolds  number  for  the 
12.4  /ym  probe.  The  plot  includes  Reynolds  numbers  computed  using  both  the 
maximum  and  minimum  probable  diameters  of  a  1 2.4  jjrr\  wire  in  accordance  with 
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the  manufacturer's  manufacturing  tolerances  (a  12.4  /jm  AWG  No.  56 
thermocouple  wire  can  vary  from  11.4  to  13.4  /ym  [Ref.  21]).  These  values 
represent  two  regimes  of  the  Reynolds  number:  0<Re<1.0  describes  a  region 
of  highly  viscous  laminar  creeping  motion,  and  the  range  1 .0< Re<  1 00  describes 
a  region  of  laminar  motion  with  strong  Reynolds  number  dependence  [Ref.  7:  p. 
309].  Table  V  summarizes  Reynolds  numbers  for  the  probe  as  a  function  of 
altitude  where  the  probe  velocity  is  4.4  m/s. 


TABLE  V.     Reynolds  Number  for  NPS  12.4  +/-  1.0  /vm  Chromel-Constantan 
Probes  from  0  to  30  km  at  v=4.4  m/s. 


ALT  (km) 

1 1 .4  //m  Probe 

12.4  /jm  Probe 

13.4  /vm  Probe 

0 

3.44 

3.75 

4.05 

10 

1.43 

1.55 

1.68 

20* 

0.32 

0.34 

0.37 

30 

0.06 

0.07 

0.07 

*  above  19  km,  where  Re  <0.4,  the  equation  for  Re  begins  to  lose  its 
accuracy;  however,  empirical  fits  to  the  data  by  Holman  [Ref.  16:p.  303] 
have  shown  that  Re  values  follow  approximately  the  same  trend  in  this 
regime. 
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REYNOLDS  NUMBER  VS.  ALTITUDE 

12.4  Micron  Probe  (V=4.4  m/s) 


4.5 


10  15  20 

ALTITUDE  (km) 


30 


Figure  18:  Reynolds  Number  vs.  Altitude  for  NPS  12.4  +/-  1.0  jjm    Chromel- 
Constantan  Thermocouple  Probe. 
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f.     Nusselt  Number 

From  Kothandaraman  and  Subramanyan  the  Nusselt  number  for 
a  cylinder  in  a  cross-flow  is  [Ref.  22:p.  91] 


Nu=C 


fu<fi 


(17) 


v   v  / 


where  C  and  n  are  empirically  developed  constants  based  upon  the  local 
Reynolds  number  where 


Re= 


<ujfi 


v   ) 


(18) 


Table  VI  contains  the  values  of  the  constants  C  and  n.  Figure  19  is  a  plot  of  Nu 
vs.  altitude  for  the  NPS  probe  based  on  the  Reynolds  numbers  shown  in  Table 
V.   Table  VII  summarizes  the  Nusselt  number  values. 
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TABLE  VI.    Constants  for  Determination  of  Nusselt 
Number  for  Eq.  (16). 


Re 

C 

n 

<  0.4 

0.891* 

0.330* 

0.4  -4.0 

0.891 

0.330 

>  4.0 

0.821 

0.385 

empirically  fit 


TABLE  VII.   Nusselt  Number  vs.  Altitude  for  NPS  12.4  +/-  1.0  /ym  Chromel- 
Constantan  Thermocouple  Probes. 


ALT  (km) 

1 1 .4  /ym  Probe 

12.4  /ym  Probe 

13.4 /ym 
Probe 

0 

1.34 

1.38 

1.41 

10 

1.00 

1.03 

1.06 

20 

0.61 

0.63 

0.64 

30 

0.35 

0.37 

0.38 
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NUSSELT  NUMBER  VS.  ALTITUDE 

12.4  Micron  Probe  (V=4.4  m/s) 


1.6 


0.2 

0 
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10  15  20 

ALTITUDE  (km) 


25 


30 


Figure  19.     Nusselt  Number  vs.  Altitude  for  NPS  12.4  +/-  1.0  /ym  Chromel- 
Constantan  Thermocouple  Probe. 
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g.    Rate  of  Convective  Heat  Flow 

Newton's  Law  of  Cooling  governs  the  rate  of  convective  heat  flow 
between  a  solid  and  its  environment  by  [Ref.  19:p.  89] 

q=hAAT,  (19) 

where 

q  =  rate  of  heat  transfer 
h  =  convective  heat  transfer  coefficient 
A  =  solid  surface  area 
AT  =  temperature  difference  between  the  solid  and  its  environment. 

h.    Convective  Heat  Transfer  Coefficient,  hc 

The  convective  heat  transfer  coefficient,  hc,  is  a  parameter  that 
describes  the  energy  transport  (both  thermodynamic  and  mass)  between  the 
probe  and  the  surrounding  atmosphere.  It  is  a  fluid  property  that  conveniently 
groups  terms  that  are  a  complex  function  of  fluid  flow  and  geometric  conditions 
of  the  probe  and  the  environment  [Ref.  19:p.  44].  It  is  generally  expressed  as  an 
empirical  form  of  the  Nusselt  number,  Nu,  by 


NuM,  (20) 

k 


54 


where  d  is  the  characteristic  length  (wire  diameter)  and  k  is  the  thermal 
conductivity  of  the  surrounding  air  in  (kcal/m-hr  K)  [Ref.  23:p.  6-3].  The 
convective  heat  transfer  coefficient,  hc,  has  dimensions  of  W/m2K. 

Previous  work  by  Brown  [Ref  12:p.  19]  and  Weitekamp  [Ref.  3:p. 
38]  has  shown  that  there  is  approximately  a  factor  of  four  decrease  in  the 
convective  heat  transfer  coefficient,  hc,  from  the  ground  to  an  altitude  of  30  km. 
The  value  of  hc  for  the  NPS  12.4  /vm  chromel-constantan  probe  exhibited  similar 
behavior;  the  actual  values  of  hc  were  different  due  to  the  different  sizes  and 
composition  of  the  respective  probes,  but  the  trend  was  the  same  as  shown  in 
Figure  20  for  the  NPS  probes  ascending  at  4.4  m/s.  Table  VIII  summarizes  these 
data. 


TABLE  VIII.    Convective  Heat  Transfer  Coefficient,  hc,  vs.  Altitude  for  NPS  12.4 
+/-  1.0 /ym  Chromel-Constantan  Thermocouple  Probes. 


ALT  (km) 

hc  for  1 1 .4  /ym 
Probe  (W/m2K) 

hc  for  12.4/ym 
Probe  (W/m2K) 

hc  for  13.4  /ym 
Probe  (W/m2K) 

0 

2973 

2808 

2665 

10 

1760 

1662 

1577 

20 

1040 

982 

932 

30 

637 

602 

571 
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HEAT  TRANSFER  COEFFICIENT,    h 

12.4  Micron  Probe  (V=4.4  m/s) 
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Figure  20.  Atmospheric  Convective  Heat  Transfer  Coefficient,  hc,  vs.  Altitude  of 
NPS  12.4  +/-  1.0  /ym  Chromel-Constantan  Thermocouple  Probe  with  Wind 
Velocity  of  4.4  m/s. 


56 


B.      TIME  CONSTANT,  r 

1.      Theoretical  Probe  Response 

The  time  for  the  thermocouple  to  complete  63.2  %  of  its  response  to 
a  step  change  in  the  temperature  of  the  air  around  it  is  the  characteristic  time,  r 
[Ref.  24:p.  566].    The  characteristic  time  is 

*=£r,  (21> 


where  p  is  the  thermocouple  density  in  (gm/cm3),  c  is  the  thermocouple  specific 
heat  capacity  in  (J/g»K),  d  is  the  diameter  of  the  thermocouple  in  (cm),  and  hc 
is  the  convective  heat  transfer  coefficient  in  (J  •  s)/(cm2  •  K)  Hence,  the 
characteristic  time  is  the  ratio  of  the  probe  thermal  storage  capacity  to  the  heat 
input  per  degree  of  temperature  difference.  [Ref.  24:p.  555] 

Eq.  (21)  shows  that  the  response  time  of  a  particular  probe  varies  only 
with  the  convective  heat  transfer  coefficient,  hc,  the  other  quantities  depend  on 
the  structure  of  the  thermocouple  and  do  not  change  appreciably  under 
experimental  conditions.  Several  factors  affect  the  accuracy  of  this  equation. 
Manufacturing  tolerances  in  wire  diameter  coupled  with  the  presence  of  a  welded 
junction  bead  that  connects  the  thermocouple  wires  near  the  center  (Figure  21) 
contribute  to  an  imprecise  prediction  of  probe  response. 
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IDEAL  WELDED 
JUNCTION  BEAD 


D  =  3xd 


4 ► 


Thermocouple  Wire 


Figure  21.    Thermocouple  Junction  Bead. 
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The  geometry  and  mass  of  this  welded  junction  bead  give  it  different 
heat  transfer  characteristics  than  those  of  the  cylindrical  wire,  and  therefore,  affect 
the  true  value  of  r.    Moffat  accounted  for  this  difference  by  the  expression 

/  n\0.375 

(22) 


t0  v  d 


where  r  is  the  true  characteristic  time,  r0  is  the  characteristic  time  when  the 
diameter  of  the  junction  bead  and  the  diameter  of  the  wire  are  the  same  (no  weld 
effect),  D  is  the  diameter  of  the  bead  (assumed  to  be  spherical),  and  d  is  the 
diameter  of  the  probe  wire.  Moffat  also  emphasized  that  to  ensure  first  order 
probe  response,  the  diameter  of  the  weld  should  not  exceed  the  wire  diameter 
by  more  than  10%.  [Ref.  24:p.  566] 

2.      Thermocouple  Characteristics 

(1)     Thermocouple  Type 

The  NPS  differential  thermosonde  used  thermocouple  rather 
than  resistance  probes  to  simplify  the  electronic  circuit.  The  smallest 
commercially  available  thermocouple  was  12.4  fjm  in  diameter,  which  was  four 
times  larger  than  those  used  by  the  Geophysics  Directorate.  The  increased 
thermal  mass  reduced  the  high  frequency  response  which  had  increasing 
importance  in  structure  function  measurements  for  small  separations  (-10  cm). 
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(2)  Material  Properties 

The  NPS  Group  used  an  ANSI  Type  E  cylindrical  wire 
thermocouple  to  make  temperature  measurements  for  all  launches  from  July  to 
November  1992.  According  to  the  specifications  of  the  manufacturer,  Omega 
Engineering,  Inc.,  this  chromel-constantan  thermocouple  was  12.4  jjm  in 
diameter.  Its  composition  was  75%  (Ni-10Cr)  and  25%  (Cu-43Ni)  which  gave  it 
a  specific  heat  capacity,  c,  of  .4367  (J/g  •  K).  Appendix  C  has  more  details  on  the 
thermocouple  composition  and  characteristics.  The  NPS  Group  chose  this 
thermocouple  since  it  was  the  smallest  available  and  because  of  its  large 
thermoelectric  response. 

(3)  Thermocouple  Geometry 

To  ensure  that  all  constants  affecting  the  time  constant  were 
valid,  we  measured  the  diameter  of  the  chromel-constantan  thermocouple  wire 
with  a  Cambridge  200,  tungsten  filament  scanning  electron  microscope  (SEM). 
The  measured  diameter  of  the  cylindrical  segment  of  the  wire  was  13.7  jjm  with 
an  unspecified  tolerance.  This  value  was  10%  higher  than  the  manufacturer's 
specifications  for  the  AWG  no.  56  thermocouple  wire.  Since  the  SEM  had  not 
been  calibrated  recently  and  the  operator's  manual  contained  no  specific 
measurement  tolerances,  we  used  the  measurement  as  a  verification  of  the 
general  dimensions  of  the  wire  and  to  estimate  the  bead  size  [Ref.  25]. 

The  "v"-shaped  wire  had  a  welded  bead  near  the  v,  which 
according  to  the   manufacturer,    should   have  been   spherically  shaped   and 
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approximately  3  times  the  wire  diameter  [Ref.  18:p.  A-31].  The  SEM  measured 
the  dimensions  of  the  roughly  cylindrical  bead  to  be  approximately  15.5  /jm  x 
15.5  jjm  x  52.3  jjm  resulting  in  a  volume  of  approximately  13,000  /vm3.  Since  the 
SEM  only  measured  in  2  dimensions,  it  did  not  provide  a  value  for  the  bead's 
height.  Photograph's  of  the  bead  showed  that  the  center  segment  of  the  bead 
appeared  to  be  cylindrical;  therefore,  we  assigned  the  bead  a  thickness  of 
15.5/ym,  the  same  value  as  the  width.  We  used  a  spherical  approximation 
because,  ideally,  this  would  be  the  shape  of  the  weld  bead  according  to 
manufacturer's  specifications  and  as  found  in  the  examination  of  numerous  other 
thermocouple  wires.  By  approximating  the  shape  of  the  bead  as  spherical  and 
using  the  standard  equation  for  the  volume  of  a  sphere, 

V=U  D\  (23) 


where  D  is  the  diameter  of  the  sphere,  the  effective  diameter  of  the  bead  became 
29  fjm.  This  diameter  was  -2.3  times  the  wire  diameter,  in  reasonable 
agreement  with  the  3:1  ratio  given  by  the  manufacturer.  We  calculated  an 
adjusted  value  for  the  bead  diameter  due  to  the  discrepancy  between  the 
manufacturer's  specifications  and  the  measured  SEM  values  for  the  wire 
diameter.  Since  the  SEM  value  was  10%  higher  than  the  manufacturer's  value 
for  diameter,  its  values  for  the  bead  dimensions  should  have  been  10%  higher 
as  well.     Factoring  this  adjustment  into  Eq.  (23),  yielded  an  adjusted  bead 
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diameter  of  26  /ym.  Using  either  the  SEM  or  manufacturer's  values,  the  ratio  of 
the  bead  diameter  to  the  wire  diameter  clearly  exceeded  the  10%  allowance  for 
first-order  effects,  therefore,  higher-order  effects  would  affect  the  probe  response. 

3.      Theoretical  Calculation  of  r 

Figure  22  shows  a  plot  of  r  vs.  altitude  for  the  NPS  12.4  /ym  diameter 
wire  thermocouple,  uncorrected  for  a  junction  bead.  The  uncorrected  r  for  the 
12.4  /ym  probe  had  a  value  of  4.2  ms  at  0  km  and  19.7  ms  at  30  km.  Figure  23 
shows  a  plot  of  t  vs.  altitude  for  the  NPS  12.4  /ym  diameter  wire  thermocouple, 
corrected  for  a  26  /ym  junction  bead.  The  time  constant  accounting  for  the  bead 
was  5.8  ms  at  0  km  and  27.0  ms  at  30  km.  The  effect  of  correcting  for  the  bead 
yielded  response  times  that  were  -30%  slower  than  for  the  calculations  with  no 
bead  correction.    Table  IX  summarizes  these  data. 


TABLE  IX.  Calculated  Thermocouple  Probe  Response  vs.  Altitude,  for  a  12.4  +/- 
1.0  /ym  Chromel-Constantan  Wire,  with  and  without  Correction  for  a  26  /ym 
Junction  Bead. 


d  (/ym) 

11.4 

11.4 

12.4 

12.4 

13.4 

13.4 

ALT 
(km) 

no 
bead 

26 /ym 
bead 

no  bead 

26 /ym 
bead 

no  bead 

26 /ym 
bead 

0 

3.7 

5.0 

4.2 

5.6 

4.8 

6.2 

10 

6.2 

8.4 

7.1 

9.4 

8.1 

10.4 

20 

10.5 

14.3 

12.1 

15.9 

13.8 

17.6 

30 

17.1 

23.3 

19.7 

26.0 

22.5 

28.8 
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TIME  CONSTANT  vs.  ALT 

12.4  Micron  Probe,    No  Bead 
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Figure  22.  Calculated  Time  Constant  for  NPS  12.4  +/- 1.0  /jm  Type  E  Chromel- 
Constantan  Thermocouple  Probe  with  no  Correction  for  a  Junction  Bead. 
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TIME  CONSTANT  vs.  ALT 

12.4  Micron  Probe  with  26  Micron  Bead 
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Figure  23.   Calculated  Time  Constant  for  IMPS  12.4  +/-  1.0 /ym  Type  E  Chromel- 
Constantan  Thermocouple  Probe  with  a  Correction  for  a  26  /ym  Junction  Bead. 
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C.       EXPERIMENTAL  VERIFICATION  OF  PROBE  RESPONSE  TIME 

1.       Empirical  Formula  For  Probe  Thermal  Response 

When  used  to  measure  local  gas  temperature,  a  bare-wire 
thermocouple  displays  physical  and  thermodynamic  behavior  similar  to  that  of 
a  cold-wire  anemometer.  A  cold-wire  anemometer  measures  rapidly  fluctuating 
flows  such  as  the  turbulent  boundary  layer  surrounding  an  object  (the  probe). 
It  consists  of  a  fine  wire  mounted  between  two  small  probes.  An  empirical 
equation  that  relates  the  wire  response  time  with  the  windspeed  is 


±=a+bjv  (24> 


where  a  and  b  are  constants,  a  being  the  y-intercept  of  the  plotted  data  and  b 
being  the  slope  of  this  data,  plotted  according  to  linear  regression  where  t  is  the 
response,  and  v  is  the  wind  velocity  [Ref.  26:p.  13]. 

2.      Wind  Tunnel  Test 

We  constructed  a  wind  tunnel  and  subjected  a  12.4  /ym  probe  (not  the 
same  probe  that  we  had  SEM  diameter  measurements  on)  to  cross-winds  varying 
in  speed  from  0  to  4.4  meters-per-second.  The  maximum  windspeed  was  4.4 
m/s  because  of  the  operating  limit  of  the  24V  DC  fan  that  was  the  wind  source. 
The  4.4  m/s  windspeed  approximated  the  average  ascent  velocity  of  a  balloon 
which  typically  ranges  from  4  to  6  m/s.   While  it  was  in  the  wind  tunnel, 
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we  illuminated  the  probe  with  an  unfocused  visible  laser  diode  to  provide  a 
periodic  instantaneous  heating  effect.  This  was  to  simulate  the  effect  of  an 
instantaneous  change  in  air  temperature  as  the  probe  moved  through  different 
air  masses  as  it  would  throughout  a  flight.  The  control  portion  of  this  experiment 
was  conducted  at  sea  level  conditions.  The  wind  tunnel  was  built  of  styrofoam 
in  order  to  minimize  the  influence  of  any  external  heat  sources  on  the 
thermocouple.  Figure  24  illustrates  the  experimental  setup.  Appendix  C  contains 
additional  details  of  the  experiment. 
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Figure  24.  Schematic  of  Wind  Tunnel  Thermocouple  Response 
Experiment  Conducted  in  Monterey,  CA  at  0  km  altitude  and  25° 
C  on  8  November  1992. 
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3.       Experimental  Results 

The  12.4  /jm  probe  had  a  measured  time  constant  of  5.6  +/-  0.2  ms 
in  a  4.4  m/s  wind  in  the  wind  tunnel  at  0  km  altitude.  Figure  25  is  a  plot  of  the 
probe  response  to  the  periodic  laser  pulse  at  this  windspeed.  It  illustrates  the 
rapid,  but  not  instantaneous,  response  of  the  probe  to  an  instantaneous 
temperature  change.  The  horizontal  distance  along  the  x-axis  between  the 
instantaneous  change  and  the  time  it  takes  the  amplitude  to  reach  the  1-1  /e  point 
is  the  time  constant  r.  Appendix  C  contains  plots  of  the  results  of  performing  this 
test  on  a  25  //m  copper-constantan  probe,  a  probe  used  frequently  in  previous 
experiments.    Table  X  summarizes  these  data. 


WIND  TUNNEL  PROBE  RESPONSE 

Type   "E"    12.4-    micron   Probe   (V=4.4-    m/s) 
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Figure  25.  Experimental  Response  of  a  12.4  jjm  Diameter  Thermocouple  Probe 
Illuminated  by  a  Periodic  Laser  Diode  in  a  4.4  m/s  Wind  in  the  NPS  Wind  Tunnel 
Experiment  of  2  November  1992.  Time  Constant  Measured  to  be  r  =  5.6  +/-  0.2 
ms. 
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TABLE  X.    Measured  Probe  Response  in  Wind  Tunnel  for  12.4  and  25  //m 
Wires. 


Wind 

(m/s) 

12.4  fjm 
Chrom-Const 

25  /ym  Copper- 
Const 

Response  Ratio 
(25//m/12.4//m) 

0 

15.7 

69 

4.4 

1.5 

7.1 

28.5 

4.0 

3.0 

6.4 

21 

3.3 

4.4 

5.6 

19.5 

3.5 

These  results  show  a  factor  of  four  increase  in  probe  response  time 
for  a  factor  of  two  increase  in  probe  diameter.  This  shows  the  effect  of  the 
thermal  mass  slowing  the  response  time.  The  quadratic  diameter  dependence 
arises  from  Eq.  (21)  because  the  hc  term  depends  on  d.  Rearranging  Eqs.  (16), 
(17),  and  (20)  yields  an  equation  of  the  form  C  x  (d2/dn)  where  C  is  a  grouping 
of  constants  and  n  is  an  empirical  constant  equal  to  0.33  for  the  Reynolds 
number  regime  of  the  12.4  //m  probe  with  a  velocity  of  4.4  m/s.  The  d2/dn  term 
then  reduces  to  d1 67  which  accounts  for  the  approximate  d2  ratio  between  the 
increase  in  the  time  constant  based  on  a  given  change  in  probe  diameter. 
Applying  this  same  rationale  to  the  spherical  junction  beads  for  the  two  different 
probe  sizes  should  account  for  most  of  the  remaining  difference. 

4.      Comparison  of  Experiment  vs.  Theoretical  Calculations 

As  previously  shown,  hc  is  a  complicated  function  of  numerous 
variables;  therefore  the  characteristic  time  constant,  r,  is  also  a  complicated 
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function  that  should  change  in  a  similar  fashion.  Table  XI  compares  experimental 

response  measurements  with  analytical  model  calculations  for  the  12.4/ym  probe. 

Appendix  C  contains  the  experimental  data  and  the  results  of  the  regression 

analysis.    We  did  not  include  the  plots  for  the  11.4  and  13.4  /ym  probes  here 

because  they  would  obscure  the  results  and  because  their  relationship  to  the 

response  of  the  12.4  /ym  probe  is  shown  in  Figures  22  and  23.    Figure  26  plots 

the  time  constant  at  0  km  as  a  function  of  windspeed.    This  shows  that  at  v  > 

1.5  m/s  windspeed,  the  analytical  model  closely  approaches  the  wind  tunnel 

experimental  measurements  for  a  12.4  /ym  wire  with  a  26  /ym  bead.    At  v  =  4.4 

m/s,  the  average  ascent  velocity,  the  curves  of  the  experimental  regression  data 

virtually  overlay  the  modeled  calculations.    At  v  <  1.5  m/s  the  plot  shows  the 

damping  effect  of  the  larger  thermal  mass.  Figure  27  shows  the  same  data  using 

Eq.  (24),  the  anemometer  response  expression. 

TABLE  XI.  Time  Constant  Comparison:  Experimental  Measurements  vs. 
Analytical  Model  for  the  NPS  12.4  +/- 1.0 /ym  Chromel-Constantan  Thermocouple 
Probe. 


V 
(m/s) 

EXPERIMENT 
Wind  Tunnel 

EXPERIMENT 
Curve  fit 

MODEL 

12.4  /ym 
No  bead 

MODEL 

12.4  /ym 
26  /ym  bead 

0.01 

15.7 

15.0 

14.7 

19.4 

1.5 

7.1 

7.5 

6.0 

7.9 

3.0 

6.4 

6.2 

4.8 

6.3 

4.4  * 

5.6 

5.6 

4.2 

5.6 

Average  ascent  velocity 
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Figure  26.  Thermocouple  Response  Times  vs.  Windspeed.  Comparison  of 
Experimental  (Wind  Tunnel)  Measurements  with  Analytical  Model  Calculations  for 
12.4  +/- 1 .0  /ym  Chromel-Constantan  Wire,  with  and  without  Correction  for  26  /ym 
Junction  Bead. 
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12.4  MICRON  CHROM-CONST  PROBE 

(l/Time  Const)  vs.  Sqrt  of  Windspeed 
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Figure  27.  1/e  Thermocouple  Probe  Response  vs.  Square  Root  of  Wind  Velocity. 
Comparison  of  Wind  Tunnel  Measurements  with  Analytical  Model  for  12.4  +/- 1.0 
jjm  Chromel-Constantan  Wire,  with  and  without  Correction  for  26  //m  Junction 
Bead. 
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For  a  12.4  //m  probe,  the  analytical  model  overestimates  the  observed 
probe  response  by  29%  when  there  is  no  correction  factor  for  the  bead.  When 
a  correction  for  the  bead,  assumed  to  be  spherical  and  with  a  diameter  of  26  /ym, 
is  made,  the  model  is  within  1%  of  the  experimental  probe  response.  It  is  not 
probable  that  the  one  correction  made  for  bead  diameter  was  sufficient  to 
improve  the  model's  accuracy  by  a  factor  of  30.  Wind  tunnel  experimental  biases 
such  as  uneven  heating  conditions,  non-uniform  windspeed,  or  other  slight 
deviations  from  standard  atmospheric  conditions  in  the  laboratory  during 
measurements  could  have  resulted  in  a  larger  difference  between  the  model  and 
experimental  results.  Coupled  with  this  is  the  uncertainty  in  the  exact  dimensions 
of  both  the  thermocouple  wire  and  junction  bead  diameters.  This  uncertainty, 
although  not  insignificant,  does  not  obscure  the  actual  decrease  in 
responsiveness,  by  a  factor  of  approximately  4  (corresponding  to  the  decrease 
in  hc  by  the  same  factor),  that  occurs  as  the  probe  ascends  from  0  to  30  km  as 
shown  in  Figures  25  and  26.  Figures  22  and  23  show  that  variations  in  the  wire 
diameter,  even  within  the  manufacturer's  specified  tolerances,  can  alter  response 
times  by  22%  under  experimental  conditions  which  equates  to  2  ms  at  10  km  and 
3  ms  at  20  km.   Table  XII  summarizes  these  data. 
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TABLE  XII.  Differences  in  Analytical 
Model  Probe  Response  due  to 
Tolerance  in  Wire  Diameters  as  a 
Function  of  Altitude  for  the  NPS  12.4  +/- 
1.0 /ym  Thermocouple  Probe. 


ALT 
(km) 

Min  t 
(ms) 

Max  t 
(ms) 

Ar 

(ms) 

0 

5.0 

6.2 

1.2 

10 

8.4 

10.4 

2.0 

20 

14.3 

17.6 

3.3 

30 

23.3 

28.8 

5.5 

D.      PROBE  RESPONSE  SUMMARY 

The  tolerance  in  the  diameter  of  the  12.4  /ym  Type  "E"  chromel-constantan 
thermocouple  wire  and  junction  bead  make  analytical  predictions  of  its  true 
response  uncertain  to  approximately  11%  (5.6  +/-  0.6  ms  at  0  km).  Additionally, 
the  3:1  ratio  of  the  bead  to  wire  diameter  results  in  higher-order  response  effects. 
The  1 2.4  //m  chromel-constantan  probe  used  by  NPS  is  four  times  faster  than  the 
25/ym  copper-constantan  probe  used  extensively  in  previous  experiments  due  to 
its  relative  mass  to  surface  area  ratio. 

The  effect  of  the  spherical  weld  junction  on  the  thermocouple  wire  was  to 
slow  the  response  time  by  30%.  Correcting  for  a  spherical  junction  bead  allowed 
the  1 2.4  /ym  probe  to  achieve  a  modeled  performance  within  1 1  %  of  experimental 
wind  tunnel  results.   A  probe  with  uniform  diameter  and  a  spherical  bead,  or  no 
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bead,  would  allow  for  even  more  precise  extrapolations  of  modeled  response 
calculations.  A  faster  probe  response  would  improve  the  vertical  layer  resolution 
as  well  as  permit  more  precise  structure  function  measurements  at  smaller  probe 
spacings.  A  probe  that  optimizes  response  by  minimizing  its  wire  diameter,  weld 
bead  diameter,  density,  and/or  its  specific  heat  capacity  would  achieve  this  result. 
A  thermocouple  made  of  the  10  jjm  wire  currently  manufactured  by 
Isabellenheutte  Heusler,  Gmbh,  in  Germany  should  be  60%  faster  than  the  12.4 
fjm  chromel-constantan  probes  used  by  IMPS,  assuming  a  proportionately  sized 
junction  bead  of  -21  jjm. 
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V.    CONCLUSIONS  AND  RECOMMENDATIONS 

A.      CONCLUSIONS 

This  thesis  investigated  two  systematic  error  effects  in  balloon  thermosonde 
systems  used  by  both  NPS  and  the  U.S.  Air  Force  Geophysics  Directorate.  It 
supports  the  existence  of  an  order  of  magnitude  solar  contamination  of  all 
thermal  probe  daytime  data  above  the  tropopause.  This  most  probably  arises 
from  a  combination  of  package  motion  (rotation  or  pendulum-like  oscillations 
producing  a  low  frequency  signal),  and  differential  solar  heating  of  the  probes 
due  to  probe  wire  oscillation.  The  high  correlation,  91%,  observed  for  the 
temperature  differences  in  the  nighttime  flight  and  the  low  correlation  obtained 
during  all  daytime  flights,  ranging  from  50%  to  82%,  particularly  near  the 
tropopause,  support  the  hypothesis  that  the  apparent  order  of  magnitude 
increase  in  optical  turbulence  is  a  solar  induced  instrument  effect. 

This  thesis  also  investigated  the  thermal  response  time  of  the  thermocouple 
used  in  the  NPS  differential  thermosonde.  The  relatively  large  size  of  the 
thermocouple  made  it  imperative  to  reconcile  its  modeled  response  as  a  function 
of  altitude  with  controlled  laboratory  tests.  The  inconsistencies  in  the  diameter 
of  the  12.4  jjm  type  "E"  chromel-constantan  thermocouple  wire  and  its  spherical 
junction  bead  limit  predictions  of  its  true  response  to  a  precision  of  11%. 
Additionally,  the  3:1  ratio  of  the  bead  to  wire  diameter  results  in  higher-order 
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response  effects.  Based  upon  calculated  changes  in  atmospheric  heat 
conductivity  the  thermal  response  time  for  the  1 2.4  /ym  chromel-constantan  probe 
should  be  approximately  9-10  ms  at  altitudes  of  10  km  and  roughly  15-17  ms  at 
20  km. 

The  geometric  characteristics  of  the  wire  and  the  junction  bead  that  govern 
probe  response  can  vary  for  each  probe.  Measuring  the  time  constant  for  a 
probe  in  a  wind  tunnel  experiment,  and  then  using  the  same  probe  in  a  launch, 
would  have  the  effect  of  calibrating  that  probe  against  the  model. 

B.      RECOMMENDATIONS 

To  eliminate  solar  effects  apparent  in  the  daytime  launches  due  to  rotation 
and  oscillation  of  the  sonde,  it  is  necessary  to  stabilize  the  sonde  as  much  as 
possible.  The  most  reliable  and  consistent  data  will  occur  when  the  probes 
maintain  a  consistent  alignment  with  respect  to  the  sun.  Possible  solutions  to 
this  may  include  the  use  of  a  gyro  mechanism  to  stabilize  the  platform,  vertical 
baffles  or  fin  stabilizers,  solar  cells  to  determine  relative  solar  angle,  and  weights 
or  anchors  to  hold  the  apparatus  steady.  To  reduce  the  effects  of  probe  wire 
vibration,  the  tungsten  wire  probes  must  be  attached  straight  across  the  probe 
supports  rather  than  in  a  flexible  configuration.  This  will  involve  a  delicate 
compromise  to  avoid  thermal  stress  and  strain  on  the  wire  and  its  contacts  and 
yet  suppress  any  vibration. 
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Comparison  of  results  from  taut  wire  probes  with  those  from  the  GP-type 
loose  wire  thermocouple  should  help  clarify  the  issue  of  probe  vibration. 
Calibrating  a  set  of  probes  in  a  wind  tunnel  experiment  prior  to  a  launch  should 
improve  consistency  in  thermosonde  temperature  measurements. 

A  faster  probe  would  improve  measurement  precision  of  the  NPS  system 
in  the  stratosphere.  A  probe  that  optimizes  response  by  minimizing  its  size, 
density,  and/or  its  specific  heat  capacity  would  achieve  this  result.  Using  a  10 
/vm  thermocouple  probe  wire,  available  in  Germany,  should  reduce  the  response 
time  by  a  factor  of  about  1.6.  Additionally,  ensuring  that  thermocouple  wire 
junction  beads  are  geometrically  uniform  would  allow  for  more  precise 
extrapolation  of  experimental  data  to  stratospheric  data.  Regardless  of  the 
probes  used,  thorough  consistency  checks  should  be  built  into  subsequent 
experiments.  Thermocouple  probes  from  the  same  lot  should  be  measured  for 
wire  diameter  and  both  bead  size  and  shape  under  a  calibrated  SEM.  These 
same  probes  should  then  be  tested  for  response  in  a  wind  tunnel  experiment, 
varying  the  laser  to  illuminate  both  the  wires  and  the  junction  beads.  Finally, 
these  probes  should  be  launched  in  a  thermosonde  package  that  has  redundant 
means  of  measurement  at  equal  distances  to  establish  an  in-flight  empirical 
database  with  which  to  resolve  atmospheric  structure  functions. 


! 
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APPENDIX  A 
DIURNAL  VARIATION  EXPERIMENT  2  DATA 

Additional  data  from  DVE  2  are  shown  below.  The  data  contained  in  this 
appendix  were  used  in  the  statistical  summary  shown  in  Table  II  and  do  not 
exhibit  any  atypical  behavior.  Figures  28  thru  33  show  the  results  of  the  4,  5,  and 
6  December  1990  flights.  Figures  34  thru  37  show  the  averaged  Cn2  profiles  for 
the  7  and  8  December  flights. 
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Figure  28.  Log  (Raw  Cn2)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  4  December  1990  (Probe  Set 
1). 
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Figure  29.  Log  (Raw  Cn2)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  4  December  1990  (Probe  Set 
2). 
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Figure  30.  Log  (Raw  Cn2)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  5  December  1990  (Probe  Set 
1). 
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igure  31 .  Log  (Raw  Cn2)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  5  December  1990  (Probe  Set 
2). 
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igure  32.  Log  (Raw  Cn2)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  6  December  1990  (Probe  Set 
1). 
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Figure  33.  Log  (Raw  Cn2)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  6  December  1990  (Probe  Set 
2). 
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Figure  34.  Averaged  Cn2  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set 
1). 
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Figure  35.  Averaged  Cn2  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set 
2). 
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Figure  36.  Averaged  Cn2  Profile  (Night),  Collected  by  the  Geophysics 
Directorate,  from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1990 
(Probe  Set  1). 
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Figure  37.  Averaged  Cn2  Profile  (Night),  Collected  by  the  Geophysics 
Directorate,  from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1990 
(Probe  Set  2). 
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APPENDIX  B 
THERMOCOUPLE  PROPERTIES 

Table  XIII  shows  the  specific  heat  capacities  of  the  materials  that 
comprise  the  type  E  chromel-constantan  probes  [Ref.  18:p.  12-107].    Chromel- 
Constantan  is  75%  chromel  and  25%  constantan.    Chromel  is  90%  nickel  and 
10%  chromium.    Constantan  is  57%  copper  and  43%  nickel  [Ref.  27:p.  133], 
Applying  these  percentages  to  the  specific  heat  capacities  in  Table  XIII  results 
in  a  specific  heat  capacity  of  C=0.4367  (J/g-K)  for  the  chromel-constantan 
thermocouple  probe. 


TABLE  XIII.     Type  E  Chromel-Constantan   Thermocouple 
Material  Specific  Heat  Capacity,  C. 


Nickel 

(Ni) 

C=0.106 

Chromium 

(Cr) 

C  =  0.110 

Copper 

(Cu) 

C  =  0.092 

Chromel 

90% 

10% 

Constantan 

43% 

57% 

Specific  heat  capacity,  C,  is  in  (J/g*K). 
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APPENDIX  C 
WIND  TUNNEL  PROBE  RESPONSE  EXPERIMENT 

Table  XIV  contains  the  wind  tunnel  measurements  on  the  12.4  jjm  probe. 

TABLE  XIV.   Wind  Tunnel  Probe  Response  Measurements  vs.  Windspeed, 
12.4  jjm  Chromel-Constantan  Probe,  8  November  1992. 


V  (m/s) 

Measured  r 
(ms) 

SQRT  of  V 

1/r 

0 

15.7 

0 

0.06 

1.5 

7.1 

1.22 

0.14 

3.0 

6.4 

1.73 

0.16 

4.4 

5.6 

2.1 

0.18 

Regression  analysis  on  this  data  yielded  values  of  a  =  0.066  and  b  =  0.054 
for  the  constants  used  in  Eq.  (24). 
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Figure  38  plots  the  response  of  a  25  /vm  copper-constantan  probe  in  the 
wind  tunnel  in  which  the  time  constant  was  measured  to  be  r  =  19.5  +/-  1.0 
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Figure  38.  Experimental  Thermocouple  Response  of  a  25  /vm  Copper- 
Constantan  Probe  Illuminated  by  a  Periodic  Laser  Diode  Measured  in  a  Wind 
Tunnel  at  v  =  4.4  m/s  in  the  NPS  Wind  Tunnel  Experiment  of  8  November  1992. 
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Table  XV  contains  the  instrument  settings  for  each  item  of  equipment 
used  in  the  wind  tunnel  probe  response  experiment. 

TABLE  XV.   Wind  Tunnel  Experiment  Instrument  Settings. 


NICOLET  PRO  30 

DIGITAL 

OSCILLOSCOPE 

HP3314A 

FUNCTION 

GENERATOR 

HP6236B 

POWER 

SUPPLY 

UKSER 
DIODE 

switch  time:  500 
/ysec 

2.02  hz 

0-24  V 

~  15  cm 
from  probe 

500  sweeps 

4.26  V  Ampl. 

120  mv 

2.02  V  Offset 
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